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Abstract 

^ i' A symmetry analysis is presented for the three-dimensional non-relativistic motion of charged 

particles in arbitrary stationary electromagnetic fields. The general form of the Lie point sym- 
metries is found along with the fields that respect them, considering non-trivial cases of physical 
interest. The restrictions placed upon the electromagnetic field yield five classes of solutions, 
expressed in terms of the vector and scalar potential. The Noether type symmetries are also 
investigated and their corresponding invariants are found. A second integral of motion, besides 

1 — 1 I 

^ , the Hamiltonian, results in three general cases. Finally, a relation between the symmetries of 

0^ ' the charged particle motion and the symmetries of the magnetic field lines is established. 

in 

^ '. 1 Introduction 
in 

o . 

■ Lie symmetry methods are one of the few systematic tools for studying nonhnear differential equa- 
tions, either ordinary or partial, revealing, if any, their integrability properties. In dynamical 
systems, continuous symmetries can be used to reduce the order of the system and in some cases 

^ ■ even completely integrate it. Along with Noether's theorem for variational problems, they may 

■ also yield conservation laws, which restrict solutions to an invariant manifold. Therefore, the di- 
rect search for symmetries in certain systems has received a great deal of attention over the past 
few decades. Another usage of symmetry analysis, which has been adopted more and more often, 
is to classify all possible symmetry groups admitted by a wide family of differential equations. 
Such results determine under which conditions a system of general form may possess one or more 
symmetries. 

Following this approach, we investigate the symmetry properties of the Newtonian motion of a 
charged particle in an arbitrary, yet stationary, electromagnetic field. Obeying the Lorentz force 
law, the system is expressed by three second order, autonomous differential equations, which in 
vector form are 

X = X X B{x) + E{x), (1) 

where B and E are smooth vector functions of the position x alone, considered in the Euclidean 
space M^, representing the magnetic and electric field, respectively, while the dot stands for deriva- 
tion with respect to time t. For B = 0, the above system reduces to a general problem of Classical 
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Mechanics, for which a complete symmetry group classification can be found in [2] in two dimen- 
sions and [3, 4, 15] in three. A symmetry analysis for the two-dimensional case of system (1), 
where the magnetic field has a constant direction in space, has also been made [7, 8], even for 
time-dependent electromagnetic fields. Therefore this case, along with trivial ones of less physical 
interest, where either one of the fields is homogeneous, will not be considered here. Our intention is 
to cover a variety of applications, in which the electromagnetic field may be quite complicated, such 
as those occuring in plasma physics and fusion devices. A series of simpler problems for particular 
choices of the functions B and E is presented in [12]. 

In the present work, we search for Lie as well as Noether point symmetries admitted by the 
above equations, leaving the electromagnetic field as arbitrary as possible, neglecting the previously 
mentioned cases. We find that the general form of Lie point symmetries is a combination of dilations, 
rotations and translations, i.e. they are linear. Starting with an unprescribed electromagnetic 
field, the symmetry condition, apart from the symmetry generator, eventually places restrictions 
on the functions B and E, too. The latter are expressed through three coupled first order linear 
partial differential equations for each field. Taking into account Maxwell's equations as well, these 
equations are solved in terms of the vector and scalar potentials, A and Thus, we determine 
the general form of electromagnetic fields of physical interest, for which systems of class (1) admit 
point symmetries, along with the general form of the symmetries themselves. More precisely we 
obtain families of potentials that are compatible with a second Lie point symmetry, besides the 
expected time translations. Since, in terms of A and the problem becomes a variational one, we 
also check which of the symmetries found are of Noether type. It turns out that, with the exception 
of one dilation, all other symmetries are variational, too. We find the correspondings invariants, 
focusing on integrals of motion, over and beyond the well-known Hamiltonian of the system. An 
important aspect is that these symmetries are also preserved by the magnetic (and the electric) 
field itself, and the induced integral for the magnetic field lines is also found. 

This paper is organised as follows. In section 2 we determine the most general form of the Lie 
point symmetry generator of equations (1), along with the constraints of the electromagnetic field. 
These are tackled in the next section, where Maxwell's equations are also imposed for consistency 
with any real physical problem. Consequently, the previous conditions on B and E are expressed in 
terms of A and and then solved. We obtain five classes of solutions, each of which corresponds to 
one specific symmetry, apart from time translations. In section 4 the general form of Noether point 
symmetries and the related electromagnetic fields are also found. The corresponding integrals of 
motion are derived in the next section. The existence of a second constant of motion, functionally 
independent of the Hamiltonian of the system, is considered, and three types of such integrals 
are obtained. Finally, section 6 introduces a comparison of the charged particle motion with the 
dynamics of the magnetic field itself, where the two systems are tested for common symmetries. 
An integral of motion for the magnetic field lines is obtained, as well. 

2 Lie point symmetries 

Using a cartesian frame of reference, system (1) can also be written as 
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where eijk is the Levi-Civita symbol. Einstein's summation convention has been adopted, assuming 
all indices from now on and throughout the rest of this paper take values from 1 to 3. 
We consider infinitesimal point transformations of the form 

t = t + e^{t,x) + 0{e^) , (3) 
Xi = Xi + eipi{t,x) + 0{e^) , (4) 



generated by the vector field 



The symmetry condition then reads 

X^^\xi - eijkijBk - E,) = 0, (6) 

whenever equations (2) hold, where X^'^^ represents the second prolongation of X. For a detailed 
description of symmetry theory in differential equations, see [1, 11, 14]. 

After substituting the second derivatives, equations (6) take the form of polynomials of third 
degree in terms of the first derivatives Xj. Being identities for all t,Xi and Xi, they finally break up 
into a larger set of partial differential equations, the so-called determining equations, 

e.,.. = 0, (7) 

"^^i^j^k ~ ^ijlix^Bi + eukCxjBi — 5ij{2^tXk — ^ImkCxiBm) - Sik{'^CtXj + f^ljmixiBm) = 0, (8) 

'^'^hxi ^ijk'PixiBk - eijkitBk - 'iixjEi - eukipi^.Bk - eijkfiBk^^ - 5ij{iu + ixiEi) = 0, (9) 

"Pitt + Vi:.^ Ej - 2CtEi - eijk^PjtBk - ^jEi^^ = 0, (10) 

where 5ij is Kronecker's delta and subscripts t or Xi denote partial differentiation with respect to 
each of them. 

Starting from the top, the first set of equations yield ^(t,a;) = fi(t)xi + w{t). As can be seen 
from the second one, mixed second order partial derivatives „ and tpi^ ^ are already considered 
equal. Imposing further integrability conditions on (8), i.e. (fi,^,.,^ = ^i.^^^^^ = "Pi.^.^.,^ we find, 
after some investigation, that fi = for non-constant B. This in turn means that 

^{t,x) = w{t), (11) 
ipi{t, x) = hij{t)xj + hi{t) (12) 

where w, hij, hi are arbitrary functions of t. 

After substitution of the above expressions, equations (9) for i = j take the form 

2/lii + (/li3 + /l3i) B2 - {hi2 +h2l)B3-W = 0, 

2A22 + {hi2 + h2i) Bs - (/i23 + hs2) Bi - id = 0, (13) 
2/133 + {h23 + h32) Bi - {his + hsi)B2-w = 0, 

from which follows the relation 

2hii = 3w + 3c, (14) 

where c is a constant. 
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The rest of (9) for i ^ j can be rearranged into two sets, by taking ij- and ji-equations in pairs. 
Adding the two equations of each pair gives the first set, 



2 [hi2 + /l21 j - (/ll3 + /l3l) Bi + (/l23 + /l32) ^2 + 2 (/in - /122) ^3 = 0, 

2 (Ai3 + A31) + (/ll2 + /l2l) Si + 2 (/133 - /ill) ^2 - (/l23 + /i32) ^3 = 0, (15) 

2 (/l23 + /132) + 2 (/l22 - /i33) Bi - (/li2 + /l2l) ^2 + (/ll3 + /^Sl) ^3 = 0. 

Before proceeding with the second one, direct inspection of systems (13) and (15) shows that 
unless 

hij = -hji, i^j (16) 

hold, then the magnetic field is at best of constant direction. This case has already been investigated 
in [8]. 

So, in light of (16), equations (15) trivially result in 

flu = /i22 = /i33- (17) 
while the remaining of equations (9), taking also into account (14), yield 

{hjkXk + hj) = (c - hjj)Bi + hijBj + Cijkhjk, (18) 

placing inevitably restrictions upon B. 

The latter show that, for a time-independent magnetic field, the functions hij, hi must be 
constants. This can be seen, if we differentiate the last three equations with respect to Xi, and then 
treat hij, hi and c as the unknowns. Hence, we arrive at an algebraic system of nine equations in 
eight unknowns (considering (16)-(17)), which is also linear and homogeneous. In the generic case 
the rank of this system is seven, and therefore we can solve for hij and hi in terms of c, which is 
an absolute constant, and the derivatives of the magnetic field. Since B does not depend explicitly 
on time, all these solutions have to be constants. 

Thus, equations (18) take the form 

dB ' 

{hjkXk + hj) = (c - hjj)Bi + hijBj, (19) 

OXj 

while (14), with regard to (17), can also be integrated, arriving at a linear expression in t for ^, 

w{t) = {2hn -c)t + ho, (20) 

where /iq is an arbitrary constant. Finally, the last set of the determining equations after substitu- 
tion yield 

dE' ( 4 \ 
(hjkXk + hj) = ylc - -hjjJEi + hijEj, (21) 

expressing the constraints on the electric field. 
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To summarize, from the relations (11), (12), (16), (17) and (20), the most general form of the 
symmetry generator is found to be 

e = (2/iii -c)t + /io, 
ipi = hiix + hi2y - h:iiz + hi, 
ip2 = -hi2X + huy + ^232: + h2, 
^3 = hsix - h23y + hnz + /13. 

We should note, however, that while all other constants are related through (19) and (21), ho is 
completely arbitrary. Thus, the symmetry algebra is spanned by dt and a second symmetry of the 
general form 

VL = {2hn - ^) 4 + + + 

which is actually a sum of dilations, rotations and translations 

, / a d d d\ d 
VL = hii 2t— + X— + y— + z— - ct— + 
\ ot ox oy oz J Ot 

, } d d\ ^ / d d\ ^ ( d d\ 

h23 z— - yTT- + ^31 X— - z— ] +hi2[ y- x— + 

\ oy oz J \ Oz Ox J \ Ox Oy J 

d d d 

hi^ + h2^ + h^ — . (24) 
Ox oy oz 

While the former, expected for every autonomous system, appears whatever the (stationary) fields 
B and E may be, the latter is admitted only when B and E satisfy equations (19) and (21). 

So, leaving the anticipated dt aside, we focus on the existence of the general symmetry vl. 
The next task to complete the symmetry analysis is to find the form of the electromagnetic field, 
respecting the symmetry condition, i.e. solve equations (19) and (21). 

3 The form of the Electromagnetic field 

The solutions to (19) and (21) describe the general form the vector functions B and E must have 
in order system (1) to admit symmetry (23). In physics, however, the electromagnetic field obeys 
certain laws, which altogether constitute another set of equations, namely Maxwell's equations. 
Gauss's and Ampere's laws simply determine the source of the electric and magnetic field respec- 
tively, i.e. the charge and current densities. But the rest two represent conditions that B and E 
have to satisfy, and are inextricably connected with every real problem (1). In the case of stationary 
fields, these are the homogeneous equations 

V • B = 0, (25) 
V X = 0. (26) 

Instead of imposing these conditions, we recall their basic consequence, that the magnetic field can 
be derived from a vector potential A = (^1,^2; ^3)) while the electric field from a scalar one 

S = V X A, (27) 
E = -V$. (28) 
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Replacing the solutions (27)-(28) in (19) and (21), provides us with the restrictions that now 
A and $ must satisfy in order to respect the Lie point symmetry condition. The equations we end 
up with, after integration, can be cast into the following form 

dAi f 2 \ df 
{hjkXk + hj) l^c- -hjjj Ai + hijAj + —, (29) 



{hjkXk + hj)— = 2(c- -hjj ]<^-k, (30) 



oxj \ 3 

where / is an arbitrary function of x and k some constant. Furthermore, the gauge invariance of 
the vector potential, A — > A + Vg, admitted by (25), can, in fact, guarantee the existence of an 
equivalent A, such that 



dAi / 2 \ 
{hjkXk + hj) ^ \~ 3^^y (31) 



for g satisfying 



cp-Vg-cg = f, (32) 

where ip is the vector with entries fi. Thus, in case of true electromagnetic fields, that certainly 
have to comply with Maxwell's equations, the solutions to (30)-(31) describe through (27) and (28) 
the ones for which system (1) has the symmetry (23). 

Consequently, we can either treat B and E only as functions entering the system, or we can 
view them as part of a bigger physical problem, that also includes Maxwell's equations. In what 
follows we focus on that second case and give the solutions A and of (30)- (31). For reasons that 
will be apparent, the potentials that the electromagnetic field comes from, besides carrying more 
information, are more convenient to use in this situation. The form of B and E can then be found 
through (27)- (28). If, however, one wishes to determine them without taking into account Maxwell's 
equations, the task would be similar. Equations (19) and (21) can be solved independently of each 
other, and obviously by the same means as (31) can. Either one of them is a system of three 
coupled first order linear partial differential equations, and in the case of the vector potential its 
explicit form is 

/ N^^l / \9Ai dAi r \ A r A 7 A 

<^i{x, y, z)— h (P2{x, y, z)— h ^z[x, y, z)—— = (c - hn) Ai + hi2A2 - /isi^s, 

ox oy oz 

, . dA2 , , dA2 , - dA2 , . , , N X 1 A , \ 

^\\x, y, z)^^ + ^2\x, y, ^)-^ + 9^3(2;, y, ^)-^ = -^12^1 + (c - /in) ^2 + /i23^3, (33) 

, , dAo, dAo^ T X A 

fiix, y,z)— h (p2{x,y,z)— h ip3{x,y,z)—— = h^iAi - /i23^2 + (0-/111)^3. 

ox oy oz 

The other two have almost the exact same form, differing only by a minor factor. Thus, we present 
in detail a way of solving (33), and then (19) and (21) can be treated accordingly. 

First of all we want to uncouple the system, which in vector form may be written as 



where Q is the 3x3 square matrix 



(cp • V)A = QA, (34) 



(35) 
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The form of the equations ahows one to do so, when Q is diagonahzable. The above matrix however 
is not, since it has two complex eigenvalues, Ai^2 = c — hn it ih, and only one real, A3 = c — /in, 
where h = (/if 2 + ^31 + ^23)^- Nevertheless, we can still separate one equation from the other 
two, by setting A = PA, where P = {u v e) is the matrix of the eigenvectors u ziz iv,e of Q 
corresponding to the eigenvalues Ai^2)A3, respectively. Then, left multiplication of equation (34) 
with the inverse of P, leads to 

i<f ■ V)A = QA, (36) 

where 

/c-hn h \ 

Q = P-^QP= -h c-hu . (37) 
^ c-hu) 

Thus, the third of equations (36) for ^43 is detached from the others. To uncouple the rest two 
we need a nonlinear transformation. Their form, as it can be seen from the matrix Q, naturally 
implies taking polar coordinates in the A1A2 plane. 

In conclusion, the following transformation in the dependent variables, i.e. the components of 
the vector potential. 



A, 



a2 = arctan f 4^ ) and A = P ^A, (38) 



as = A3 

where P is defined through the eigenvectors of Q, uncouples the three equations of (33). Indeed, 
the latter take the form 

9ai \9ai dai 

<^i{x, y, z)^— + (P2{x, y, z)^— + (fsix, y, z)^:— = (c - /ill) ai, 
ox oy oz 

y, z)^- + 'P2{x, y, z)^- + ip3{x, y, z)^- = -h, (39) 
ox oy oz 

<^i{x, y, z)-Q^ + 'P2{x, y, ^)-^ + ^3{x, y, z)-q^ = (c - /in) 03- 

Now, each equation can be solved independently by the method of characteristics. Even more 
conveniently, all three of them are essentially (the first and the third exactly) the same. 

This means that three out of four characteristic equations for the above equations are common 
in each case, forming, in light of (22), a linear dynamical system, 

=^ = Hx + h, (40) 
as 

where H is the 3x3 square matrix (hij) and h the column vector {hi, /12, /13). Since the matrices 
Q and H share the same eigenvectors, the homogeneous counterpart of the above equations can 
be easily resolved, very similarly to (34). Only now, due to the term hi, we need to make a slight 
adjustment. More specifically, by introducing the new independent variables, 

x = -\/a;2 + i/2 

y = arctan^^^ and x = {x — k) , (41) 
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where P is the same as above and k = {ki, ^3) is a constant vector soon to be defined case by 
case, we put system (39) into a much simpler form, 

'^i(^)^ + '^2(y)^ + ^3(2)^ = (c - hii) ai, 

<^i(x)^ + My)-g^ + ^3(2)^ = -h, (42) 

^l{x)—- + (P2{y)^:^ + V3{Z)^^ = (C - /ill) 03. 

ax ay az 

The solution of this system can be found easily. Notice that x, y, z describe cylindrical coordinates 
in the xyz-space, which in turn is a linear transformation of the original space. 

Summing up, transformation (38) uncouples the equations of the system, the same way (41) 
yields its characteristics. Under these two changes of variables we arrive at system (42) and then 
solve it. After the solution is found, the inverse transformations of (38), (41) give the solution of 
the original system (33). Though this treatment is general, we distinguish four characteristic cases, 
depending on the form of the matrix P and the vector k, as well as a last and more trivial one, 
where no transformation is necessary at all. 

3.1 The case hu 7^ and /i23 or /131 7^ 

In this general case the eigenvectors of Q can define the matrix 

/ /ll2/i23 -hh^i /l23 (/ill +/ii3)'^ 

huhsi /1/123 ^31(^31 + ^23)^ ' 

\_^hl + hl,) hu{hl + hj,)y 



P 



(43) 



which for this particular choice is orthogonal, P ^ = P^ , while k = —H ^h, making (40) homo- 
geneous in terms of x,y,z. So, the new independent variables, according to (41), are 



X = y/{x - ki)"^ + {y - k2Y + {z - h^Y - P 

y = -ctan { ^l^'^ - - - , , ^ (44) 

\h\2h2?, {x - ki) + huhsi [y - ^2) - [h^i + h^^) [z - k^) j 

2 = [/i23 (2; - k\) + /i3i (2/ - k2) + h\2 {z - ks)] 

Together with (38), they transform system (33) into (42) for ^i(x) = hiix,(p2{y) = —h,ip^{z) = 
/ill J. Its solution is now easily found to be 

ai = z'^ii Fi(ci,C2), a2 = y + F2 (ci,C2) , 03 = z'^n ^3(01,02), (45) 

where Fi,F2, F3 are arbitrary functions of the characteristics 

z 

ci = — , 

x' (46) 
C2 = hiiy + hlnz. 
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Taking the inverse transformation of (38), the previous solution can be expressed back in the 
original components of the vector potential, arriving at 

A = z^~^PR{y)F{ci,C2), (47) 

where F = (Fi, F2, -F3) is an arbitrary vector function of ci and C2, coming from F = (Fi, -^3)) 
and R is the rotation matrix around the third axis. 



Ri-d) 



^ cos 1? — sin 0^ 

sini? cosi? I . (48) 
V 1^ 



In terms of the original coordinates, solution (47) is further reduced, though some of the transformed 
variables are still kept for simplicity, 

Ai = z''n {[hix - ki) - h23z\Fi (01,02) + [/isi {z - h) - hu [y - k2)] F2 {01,02) + /123Z-F3 (ci,C2)} 

_ c 2 ^ 

A2 = zi^ii {[h {y - k2) - hsiz] Fi {01,02) + [/ii2 {x - ki) - /i23 {z - ks)] F2 (ci,C2) + hsizFs {01,02)} 

^3 = z'^ii {[h {z - k-i) - hi2z\ Fi {01,02) + [/i23 {y - k2) - hsi {x - ki)] F2 (ci,C2) + h^zF^ {01,02)} 

(49) 

Equation (30) on the other hand is very similar to either one of system (42), and can be solved 
likewise. So, the scalar potential is 

^^^A^~')G{o^,02), c^/in ^50^ 
^y + G{oi,02), o = hii 

up to some additive constant, where G is arbitrary. 

Since ki^s are quite complicated expressions of /ij's, it is preferable in this case for practical 
purposes to consider the latter in terms of the former, and thus, to conclude that the potentials of 
the form (49)-(50) yield symmetry (23) for hi = —hijkj. 

3.2 The case hu = and /i23 or /131 ^ 

Here we may use again the matrix P given in (43), but, since H is no longer invertible, we define 
k = ^Hh. The new coordinates resemble the ones of the previous case, 

x = y/{x - fei)2 + {y- k2)'^ + {z- ks)^ - P 

y = arctan ( ^[^23 {y - h) - {x - h)] ^\ ^^^^ 

\h12h23 {x - ki) + /ii2/i3i {y - k2) - (/131 + /123) {z - h) J 

i^23X + h-iiy + hi2z) 

where ki = {hi2h2 - h^ih^) /h? ,k2 = {h23h3 - h^hi) / h"^ , k^ = {h3ihi - h23h2) / h? , but now 
^i(x) = Q,^2{y) = -h,ip3{z) = h3, where /13 = (/i23/ii + /i3i/i2 + /ii2/i3) /h. So, 

ai = e~T^yFi{o3,Ci) , 02 = ^ + ^2(03,04), as = e~'^^F3 (03, C4) , (52) 
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where the arbitrary functions Fi,F2, depend now on the characteristics 

^ (53) 
C4 = + hz, 

that actuahy being the essential difference with the first case. 
Thus, the vector potential in the original reference frame is 

A = e-^ypR{y)Fic3,Ci), (54) 

or more simply 

Ai = e~hy {[h{x - ki) - h23z\ Fi (03,04) + [/131 {z - k^) - hi2 {y - /C2)] -F2 (£3,04) + /i23-^3 (c3,C4)} 

A2 = e~T^y {[h{y- k2) - /i3i2] Fi (03, C4) + [hi2 {x - ki) - h23 (z - ks)] F2 (03, C4) + h^iFs (03, C4)} 

^3 = e~hy {[h{z - kz) - hi2z\ Fi (03,04) + [/i23 {y - ^2) - /131 {x - ki)] F2 (03,04) + /ii2i^3 (£3,04)} 

(55) 



and the scalar one. 



2c ~ 



, e--s^G(03,04), 0/0 
$ = <( ^ ^ ^ (56) 

^y + G(c3,04), = 
where F and G are now arbitrary functions of 03 and 04. 

3.3 The case hn ^ and /123 = /i3i = 

Now the third equation of system (33) is already isolated and so it can be solved independently, 
yielding Aj, directly. This means that simply P = I, which is actually the only difference with the 
first case. All the results obtained there can be reproduced here by simply making this substitution. 
In particular, setting again k = —H~^h, the form of system (42) in the new variables 



X = \l {x- kif + {y - k2) 
y = arctan 



y^] (57) 

\x — ki 



Z = Z — k 



3 

and consequently its solution remain the same. Therefrom, we end up with the following expression 
for the vector potential 

A = J^-'i2(y)F(oi,C2), (58) 

where the arbitrary function F depends again on oi and 02 , but for the variables x, y , z defined 
above, noting that h = hi2 in this case. Thus, in terms of the original coordinates the above 
solution is written as 



(59) 





= {z 


- ks)"^^ 


—2 

[{x - ki)Fi (01,02) 


- (y - A;2)-F2 (ci,02 


A2 


= {z 


c 

— ks) '^11 


[{y - k2)Fi (01,02) 


+ {x- ki)F2 (01,02 


A3 


= {z 


- fes)^ 


F3 (C1,C2) 
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where 

z- k3 
\J{x- kif + {y - k2Y 



C2 = /ill arctanl — ) + hi2 In iz — k^) , 

\x- ki' 



(60) 



and ki = (/ii2/i2 - huhi) / (/ifj^ + /if2) ,k2 = - {huhi + /iii/i2) / {hfi + h\2) , ks = -h^/hu. The 
scalar potential accordingly is 

|(^_A:3)K^-0g(ci,C2), Cy^/ill ^g^^ 

[-^ln(z-/c3) + G(ci,C2), c = /iii 

where G is again arbitrary in its arguments. 

When hi2 = 0, then the characteristic C2 is just y. Thus, solution (58) simply reduces to 
A = [z — /c3)''ii F{ci,y). In other words there is no need for transformation (38) or (41) at all, 
because, as we can see from (33), the system is completely uncoupled. 

3.4 The case hu = /i23 = /j-si = 0, hu 7^ 

This is actually a combination of the last two cases. The new variables are now defined for P = I, 
as in the previous case, but k = -^Hh as in the second case, meaning 



X = \l {x- kif + {y- /C2) 
y = arctan 



( (62) 

\x-kij 
"z = z 

where ki = h2/hi2,k2 = — /11//112 and k^ = 0. The form of ipi in terms of the above variables is 
also the same as is in the second case for h = hi2 and h-^ = h^. Thus, the vector potential is 

A = e"^^i?(y)F(c3,C4), (63) 

or 

Ai = e '^12 ^"-'^1^ [(2;-/ci)Fi(c3,C4)-(y-fc2)F2(c3,C4)] 

A2 = e"^ ^^M^) [(y - k2) Fi (C3, C4) + (X - fei) F2 (C3, C4)] (64) 

A3 = e '•12 ^^"''I'^-Fs (c3,C4) 



- arctanl - — ^ 



and the scalar one, 

^^|e-^--™G(c3,C4), c/0 ^^^^ 

lhT;arctan(fE|) +G(C3,C4), c = 

where functions -F and G are arbitrary in the indicated arguments. The latter are defined again 
through (53), but for the variables introduced in this case, meaning 

C3 = \J{x-kif + {y-k2f, 

. . (y-^A^i. ^^^^ 

C4 = /i3 arctan — + ni2^;. 

\x-kij 
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3.5 The case hij = 0, /ii ^ 

Now no transformation is needed at all. Equations (33) are uncoupled, and each one of them can 
be solved independently. The solution for the vector potential is 

— X 

A = e^i F{h2X - hiy,hsx - hiz) , (67) 

while for the scalar 

e^i Gihox — hiy,h%x — hiz) , c 7^ , , 

^ ^ I V ^ itf, .5 1 ; , T- ^gg^ 

— ^x + G{h2X — hiy, h^x — hiz) , c = 
for some arbitrary vector function F and a scalar one G. 

Remark 3.1. In all of the above cases, the constants hij, hi (or ki in the first case) and c, that 
appear in the symmetry vl, are also present in the expressions for A and In this way they 
define both the electromagnetic field and at the same time the symmetry admitted by system (1). 
In other words, the cases presented in this section represent families of potentials, which, apart 
from time translations, are compatible with a second Lie point symmetry of the form (23). So, vl 
is actually one symmetry, and along with dt they generate two-dimensional symmetry algebras for 
each case respectively. Further restrictions on the arbitrariness of A and <I> can yield symmetry 
algebras of higher dimensions, but these cases will not be considered here. 



4 Noether point symmetries 

In retrospect, the introduction of the potentials of the electromagnetic field leads immediately to a 
different viewpoint regarding equations (1): They form an Euler-Lagrange system, stemming from 
the function 

L{x,x) = ^x^ + x ■ A{x) -<^{x). (69) 

Subsequently, the question of Noether symmetries, providing us with first integrals of motion, 
naturally arises. Unlike the case of Lie point ones, now Maxwell's equations are obviously implied 
from the beginning. 

Considering that Noether point symmetries are a subset of the Lie point ones, already found, 
we don't have to employ a new query from the beginning. We may check instead if and when the 
general form of the symmetry generator (22) is of Noether type. Following [13], we note that the 
vector field (5) defines a Noether point symmetry, if a function f{t, x) exists, such that 

Substituting the time-independent Lagrangian given in (69) and the components of the gener- 
ator from (22), results in a second degree polynomial in terms of the velocity components. This 
equation must hold identically for all t,Xi,Xi and so every coefficient of this polynomial must vanish, 
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yielding the following five equations 



c = 0, (71) 

dAi 2 df 
{hjkXk + ^i) ^ + 3 ^jj^i - ^ij^i = C^2) 

2 dj 

The first equation clearly rules the dilation tdt out as a noether symmetry candidate. On the other 
hand, the integrability conditions of (72) and (73), fx^xj = fxjXi and ftxi = fxit-, which guarantee 
the existence of /, expressed in terms of B and lead back to (19) and (21) again for c = 0. 
Furthermore, from (72) and (73), we can easily deduce that f{t,x) = kt + f{x), where A; is a 
constant and / arbitrary. Thus, we recover the restrictions (29) and (30), for c = 0, which, going 
the other way round, were derived, when (19) and (21) were integrated. 

So, the Noether symmetry condition, apart from (71), does not involve any other constraint. Of 
course the very expression of the Lagrange function, requires Maxwell's equations from the start. 
For a real physical problem these would be implemented in the Lie symmetry well. On 

this ground, the conditions placed upon the electromagnetic field, in terms of A and are the 
previous ones for c = 0. Once again the gauge invariance of the vector potential can be used as in 
section 3, where the restriction (32) for c = now reduces to (p ■ = f. In conclusion, equations 
(1) admit Noether point symmetries of the general form 

^ = 2hnt + ho, 

ifi = hiix + hi2y - h^iz + hi, ^^^^ 
^2 = -hi2X + hiiy + h23Z + /12, 
^3 = h^ix - h23y + hiiz + /l3 

for electromagnetic fields coming from the potentials described earlier in the previous section, by 
setting c = 0. 

We emphasize again that only /iq is completely arbitrary, corresponding to the symmetry dt- 
The rest of the constants define a second Noether symmetry, 

9 d d d 

VN =2/iiit— + (^1— + 932Tr + 953^, (75) 
at ox oy oz 

which the system has for a particular form of the electromagnetic field, specified again by them. 
While the invariant that corresponds to the first symmetry dt is the well known Hamiltonian of the 
system, a second integral of motion may arise from vn- 

5 Integrals of motion 

According to [13], the integral of motion, which via Noether's theorem corresponds to a symmetry 
of (69), is given by the following equation, 

dL ~ 

I = ^L+{ifi- ixi) Q-,-f- (76) 
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Substitution of the Lagrangian of the system yields 

/ = ■ + A) - e Q i2 ^ - /. (77) 

Equations (72) and (73) for ipi = and ^ = 1 triviahy result in a constant function /. Thus, 
from the symmetry dt, we recover through (77) the integral 

H{x,x) = ^x^ + <!>{x), (78) 

which is commonly used as the Hamiltonian function of the system, expressing the particle's energy. 

Before finding the integral that corresponds to the second symmetry (75), we comment that a 
second constant of motion would be of real value, if it is functionally independent of the already 
known Hamiltonian and in involution with it with respect to the corresponding Poisson bracket. 
For this to be the case we set hn = 0, and also require / to be time- independent, that is /c = 0. 
Then, using the gauge transformation A — > A + \/g of section 3, / = / will not enter at all in 
(77), since equation (32) for c = yields (p ■ Vg = f. Thus, we restrict our attention to a linear 
integral, 

I{x,x) = (f-{x + A), (79) 

either in terms of the velocities or the canonical momentums pi = Lj,. = Xi + Ai. This is actually 
a generalisation of the relevant two-dimensional results in [5] and [9], and not included in Lewis's 
search for quadratic invariants [10]. It is worth noticing that in this way the scalar potential defines 
the form of the Hamiltonian, while the vector potential the form of the second integral. 

We distinguish the following three cases, which correspond to the second, fourth and fifth one 
described in section 3 for c = A; = 0. 

5.1 The case /123 or h-^ ^ 

The vector and scalar potential in this case are given in (55) and (56), respectively, for c = A; = 0. 
The corresponding integral is 

i = ihviy - h-siz + hi)x + {-hi2X + h23Z + /12) y + {hsix - /i232/ + ^3) z 

- h?x^F2 {x, hsy - hz) + /i/isFa (x, /isy - hz) (80) 

where x,y,z are given in (51). 

5.2 The case /i23 = /iai = 

The vector and scalar potential in this case are given in (64) and (65), respectively, for c = k = 0. 
The corresponding integral is 

I = {hi2y + hi)x + {-hi2X + h2)y + h^z - hi2X^F2{x, h^y - h^z) + hsF-six, h^y - /112F) (81) 

where x,y,z are given in (62). 
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5.3 The case hij = 0, /ii ^ 

The vector and scalar potential in this case are given in (67) and (68), respectively, for c = k = 0. 
The corresponding integral is 

/ = hix + /i2y + h^z + hiFi{h2X — hiy, h^x — hiz) + h2F2{h2X — hiy, h^x — hiz) 

+ h3F3{h2X - hiy,h3X - hiz) (82) 

Remark 5.1. In all of the above cases the scalar potential reduces to an arbitrary function of the 
related characteristics. For each case there also exists a suitable coordinate system, where the 
vector potential is an arbitrary function of the corresponding characteristics, as well. 

There is another justification for focusing on these types of integrals, which can be apparent 
when investigating the dynamics of the magnetic field itself. The latter offer a better insight into 
the magnetic fields found in section 3. 



6 Magnetic field lines 

In this section we compare, in terms of symmetries, system (1), which describes the particle's orbit, 
with the system of equations 

f = Bix), (83) 

which describes the integral curves of the magnetic field B, commonly known as magnetic field 
lines. The independent variable here, denoted by r, is related to the line element of these curves. 
From the physical point of view, since t and r carry very different meanings, such a comparison 
can only be made on the ground of spatial symmetries independent of them, i.e. generated by a 
vector field of the form X = ipi{x)dxi- 

For system (1) these can be recovered from section 2 for ^ = 0, meaning c = 2/iii and /iq = 0. 
So, in this case, the coefficients ipi are restricted to the form (12) together with (16) and (17), while 
the conditions (19) for the magnetic field are 

On the other hand, by prolonging X up now to the first order derivatives x-(r), the symmetry 
condition for system (83) 

(x- -B,i) =0 (85) 

on its solutions, leads simply to 

dBi dipi 

without predefining the form of the symmetries. 

From the above we conclude that, for a given magnetic field, which satisfies equations (84), the 
symmetries of (1) cannot meet the requirements (86) for system (83) when hn ^ 0. All the same, 
for /ill = 0, where conditions (84) and (86) are identical, system (83) may have more symmetries 
than the linear ones of system (1). Thus, in general, not every symmetry of the charged particle 
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motion is a symmetry of the magnetic field lines and vice versa. If, however, we limit our choices 
to symmetries of the form 

d d d 

V = {hi2y - h^iz + hi) — + {-hi2X + h23Z + — + {h^ix - h23y + /13) (^7) 

ox ay oz 

then these are preserved by both systems for the same magnetic field B. 

Actually, we have already encountered v: it is a Noether symmetry of (1), that corresponds 
to the integral (79) described in the previous section. However, (83) is an Euler-Lagrange system, 
too, for the Lagrangian function 

L„^{x,x') = x' ■ A{x). (88) 
The Noether symmetry condition for these equations, 

dLm ^ difi dim _ df 
* dxi dr dx[ dr ' 

is also satisfied by v, giving = and the previous constrains on the vector potential, i.e. equations 
(29) for c = /ill = 0. Using the same gauge transformation to determine A, the corresponding 
integral, ^PiL^i — f, for system (83) becomes 

/m = V • A. (90) 

This is a projection of the integral (79) on the original configuration space M^, and therefore, 
according to the form of the symmetry and the vector potential, it, too, is seperated into three 
cases, which are the ones of the previous section without the velocities. 

In conclusion, whenever the motion of the charged particle is confined in the hypersurface 
cp ■ {x + A) = const., the magnetic field lines lie on the surface cp ■ A = const.. In this case both 
systems enjoy a symmetry of the form (87). 

Remark 6.1. The system of the magnetic field lines is integrable, when it admits the above symmetry 
V, as shown, for example, in [6]. 



7 Conclusions 

We have found five classes of stationary electromagnetic fields in terms of the potentials, which yield 
a second Lie point symmetry, besides time translations, for the three-dimensional non-relativistic 
charged particle motion. The analysis showed that in the case of non-homogeneous arbitrary fields 
the only possible symmetries are linear, consisting of dilations, rotations and translations. Besides 
time dilations, they have all proven to be of Noether type, too, without any further restrictions on 
the potentials of the electromagnetic field. The corresponding invariants have been constructed, in 
particular focusing on three cases, where the integrals are functionally independent of the Hamil- 
tonian and in involution with it. Thus, a total reduction of four can be made, and then further 
investigate the system in only two variables. 

Another aspect of these three cases is that the same symmetry is also admitted by the system 
of the magnetic field lines, yielding an integral of motion for the latter, too. Time-independent 
magnetic fields lying on a surface, with some geometrical (usually axial or helical) symmetry are 
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very often used to describe the equilibrium state of plasma configurations in the context of ideal 
Magnetohydrodynamics. The results obtained in this work may be useful for the determination 
of such symmetric magnetic surfaces, by finding new solutions of the Grad-Shafranov equation, 
for example. In addition, the comparison between the two integrals, (79) and (90), could also 
relate the behaviour, and possibly the confinement, of the charged particle to the dynamics of the 
magnetic field. This relation could be further analysed, after reduction of the particle's trajectories 
and integration of the magnetic field lines. 

A desirable development of the present work would be the treatment of the non-autonomous 
case, where there are no known integrals in general. Another possible extension is the investigation 
of complete integrability of system (1), i.e. considering subcases of those presented in section 3, 
admitting more than one symmetry. 
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